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SUMMARY 

Coupled theoretic sheaths given in this report approximate thermionic diodes. The 

sheaths result from particle emissions of electrodes affected by sheath fields and of 

cesium plasmas ionized by thermal electrons. By matching their net current densities, 

the analysis pairs collisionless emitter and collector sheaths to simulate diodes. Midway 

between the electrodes in both space and temperatures lie the plasmas at 1700° to 2000° K 
12 13 

with 10 and 10 electrons per cubic centimeter. The emitters and collectors differ by 
400° to 800° K. For these conditions, the computed outputs for the converters range 
from 0. 01 to 1 watt per square centimeter. Although such powers fall below desired 
performance levels, they provide low-current checks for thermionic -diode experiments. 

Sheath properties for the present systems and for electrodes in plasmas with raised 
electron temperatures agree with equilibrium results. But the theory applies only for 
equilibria, near -equilibria, and relatively low net transports. Overall sheath character- 
istics correlate with electrode and plasma conditions. One helpful parameter is a hybrid, 
the emission Debye length, » 6. 9 (Tg/N e p)'*' // ^. This variable correlates data for 
these sheaths better than the plasma Debye length. Effective sheath widths for the present 
study lie between 1. 7 and 2. 6 emission Debye lengths. 

The independent variables for the present graphic correlations of sheath properties 
involve only electrode and plasma parameters. Because one states these conditions at the 
outset, the graphs allow immediate estimates of sheath characteristics - without the 
complex iterative calculations. 


INTRODUCTION 

This study estimates several potential profiles that might develop in a cesium therm- 
ionic diode. The theoretic results detail pairs of collisionless emitter and collector 
sheaths (ref. 1) but leave the plasma losses to other analyses (refs. 2 to 4). In the model, 
the sheath separates a thermally ionized plasma (ref. 5) from an emitting electrode 
(refs. 6 to 9). Sheath fields alter the emission barriers (refs. 10 and 11) or effective 
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work functions of the electrodes. This field phenomenon requires a special virtual 
electrode (ref. 1) that allows both positive and negative particles to undergo the same 
absolute potential drop through the sheath. For usual thermionic -diode work, the ac- 
celerated particle from the electrode changes less in absolute potential energy than the 
retarded particle. If both charges undergo the same absolute potential change in the old 
sheath model, net currents occur in calculations for equilibria; of course, such occur- 
rences violate the second law of thermodynamics. Because of its self-consistent sheath 
drops, the present theory allows no net currents in isothermal systems. But when the 
electrode temperature rises above that of the plasma, net negative charge flows from the 
electrode; then it is an emitter. If the electrode temperature falls below that of the 
plasma, the electrode is a collector. 

By setting the plasma temperatures between those of the electrodes in a system and 
by matching the entering and leaving current densities, the present calculation method 
yields pairs of emitter and collector sheaths. Although the analysis ignores energy 
transfer from emitter to collector, placing the plasma temperatures midway between those 
of the emitter and collector reduces the differences in heat transfer. As previously 
stated, the theory applies only for relatively small net transports. Thus, with electrode 
and plasma temperatures and electron concentration specified, this method yields work 
functions and sheath data for well coupled emitters and collectors. 

Such detailed internal data for diodes are needed. Analyses of outputs alone fail to 
reveal mechanisms for transport and ionization in thermionic converters (ref. 2). Theory 
and experiment must probe the interior of the diode to provide design criteria. Identities 
of controlling processes emerge when interelectrode measurements compare with theoretic 
conditions, like those presented in this report. These processes limit diode currents and 
voltages in the plasma and at the electrodes. Reducing these output restrictions depends 
on knowing and using the internal functions of the thermionic converter. 

One example of this approach involves the ” isothermal-diode theory” (ref. 12). That 
model provides equilibrium data to compare with diode operations and diagnostics like the 
emission spectroscopy of the plasma gap (ref. 13). After an isothermal check and 
calibration of the diode and sensors, results for power production have firm bases and 
more meaning. In addition, internal measurements track the changes of mechanisms as 
the diode departs from equilibrium. But the sheath model used in the isothermal-diode 
theory merely adds simple space-charge curves (ref. 14) for ions and electrons. The 
present sheath model applies more accurately to diode equilibria. And the near- 
equilibrium and small-transport interpretations (ref. 1) allow this theory to approximate 
some heterothermal diode systems. 

Examples given in this report combine emitter and collector sheaths with like current 
densities to form diodes. The electrodes bracket cesium plasmas with several sets of 
conditions presented in table I. 
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TABLE I. - SHEATH PARAMETERS 


Cesium plasma 

Electrodes 

Electron 

concentration, 

cm -2 

Temperature, 

°K 

Emitter 

temperature, 

°K 

Collector 

temperature, 

°K 

10 13 

2000 

2400 

1600 

CO 

i-H 

o 

i-H 

2000 

2200 

1800 

10 12 

1800 

2000 

1600 

10 12 

1700 

2000 

1400 


Sheath calculations for work functions from 1.5 to 5.0 volts provided bases for 
interpolations to pair emitters and collectors. These graded solutions also yielded 
properties to compare with the general correlations of reference 1. Thus, overall 
sheath characteristics for the systems of table I appear with data for equilibria (ref. 1) 
and for electrodes in plasmas with elevated electron temperatures (ref. 1). After this 
check of the results, interpolations revealed possible combinations of electrodes for 
thermionic converters. Work functions, net current densities, sheath drops and widths, 
and potential, distance diagrams indicate theoretic conditions in the selected cesium 
diodes. 

In addition to this report and the theoretic background given in reference 1, the 
FORTRAN statement of the computing method for this sheath model is presented in 
reference 15. 


THEORY 

Reference 1 gives details of the model for emission sheaths, but to complete the 
description of potential profiles in cesium diodes, this report presents the theory also. 
However, in this description the equations appear in simple unexpanded forms. 


Assumptions 

Sheaths for emitting electrodes submit readily to analyses based on existing theories 
as the result of three assumptions: First, confining the study to plane sheaths greatly 
simplifies calculations. Second, using equations for equilibrium phenomena allows 
concise mathematical statements of emission (refs. 6 to 9) and plasma (ref. 5) processes 
In the present model, the contact ionization probability accounts for ion generation at the 
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electrodes (ref. 9). The contact ions and the electrons escape the electrode surfaces by 
thermal emission (refs. 6 to 9) into the sheath. Near the electrodes, though, sheath 
fields alter the emission barriers (refs. 10 and 11). The resulting effective work 
functions yield emitted current densities different from those for emissions with no ex- 
ternal fields. Across the sheath from the electrode, plasma atoms ionize by colliding 
with thermal electrons (ref. 5). But the third assumption, ruling out collisions in the 
sheath, removes great complications. Without collisions, all charge transport involves 
only acceleration or retardation in the sheath fields. Of course, the atoms move back and 
forth between the electrode and plasma unaffected by the sheaths. 

With the component theories identified by the names of their authors, the description 
of the sheath model is clear and brief: This study analyzes plane collisionless electron 
and positive -ion sheaths between Saha plasmas and Schottky versions of Richardson- 
Dushman, Saha-Langmuir electrodes. Although the contributing theories derived from 
restrictive models, experiments reveal that they perform well in many nonequilibrium 
situations (refs. 2, 11, and 16). These findings encourage the use of the present sheath 
model for electrode, plasma systems near equilibrium or with relatively small net 
transport (ref. 1). This might be termed the fourth assumption. 


Parameters 

A sheath solution usually begins with a listing of the controllable or independent 
variables. In a simple system, these parameters describe the plasma and the electrode, 
not the sheath. 

Plasma: 

Electron and ion concentrations (equal) 

Electron, ion, and atom temperatures 
Ionization potential 
Vapor pressure of plasma chemical 
Minimum mean free path 

The last two properties enter computed tests to prevent continuous removal of the plasma 
chemical by condensation on electrodes and to assure essentially collisionless sheaths. 
Electrode: 

Temperature 
Work function 

In some systems, the work function depends strongly on the arrival rate of the plasma 
chemical at the solid surface and the electrode temperature as well as its composition. 

In this report, the work function serves as a parameter, but it might be specified in 
accordance with plasma, sheath, and electrode properties. 
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Saha-Langmuir Null Point 


With the independent variables and assumptions given, simple calculations predict 
an isothermal electrode, plasma system with no sheath and, of course, no net current. 
In this state, each particle flux from the electrode equals the random current density of 
that particle in the plasma. The plasma results from ionization by thermal electrons 
(ref. 5): 
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(The symbols are defined in appendix A. ) The plasma indicated in equation (1) equili- 
brates with the emitting electrode at the Saha-Langmuir null point to produce a state 
described by equations (2a) to (7a). 
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Equation (2a) comes from references 6 to 8, equations (3a) and (4a) derive from refer- 
ence 9, and equations (5a) to (7a) depend on reference 5. 

The Saha-Langmuir null point defines the plasma electron potential as that at the top 
of the electrode emission barrier. In this report, particle potential means motive 
(refs. 1 and 17). Thus, as figure 1 shows in a motive diagram, the electron potential in 
the plasma stands above that in metal by a voltage equal to the work function. 
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When the work function falls below that of the Saha-Langmuir null point, excessive 
electrons escape the solid, and an electron sheath separates the electrode from the 
plasma. If the work function climbs above that of the Saha-Langmuir null point, the 
metal emits too many ions, and a positive -ion sheath intervenes. In either event, 
though, the sheath arises to prevent net current in the isothermal system. This balance 
maintains because the sheath accelerates one emitted charge and retards the other. Of 
course, this process reverses for particles entering the sheath from the plasma. 

A sheath imposes an electrostatic field on the electrode; the field changes the 
emission barrier (ref. 10). Figure 2 superposes an external field on the image potential 
in the usual manner. By this mechanism, the emission barrier or effective work function 
decreases for the particle accelerated from the electrode. 

At the Saha-Langmuir null point, though, the effective and normal work functions and 
the plasma potential are all the same voltage above the electrode Fermi level: 


, xT /3 , m , „ T 3, 27rm k 
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( 8 ) 


Of course, changing any of the independent variables in equations (2) to (7) destroys the 
Saha-Langmuir null point, and immediately, the relation of the electrode to the plasma 
grows complicated. To simplify the problems and eliminate errors, this analysis uses 
an effective electrode introduced in reference 1. 


Virtual Schottky Electrode 

The virtual Schottky electrode (ref. 1) is a theoretic construct that replaces the over- 
all functions of the electrode surface in the present sheath model. Because the potentials 
and dynamics of the electrode surface evade definition, such substitutions of gross effects 
for detailed local mechanisms are necessary as well as convenient. For example, the 
traditional thermionic emitter (ref. 12) looses one charge at the normal work function and 
the other at an effective work function lowered by the sheath field. Of course, these two 
potential levels occur neither on the electrode surface nor at the same place in space. 

But both work functions act as properties of the solid surface in the model. Furthermore, 
the traditional thermionic emitter uses the image potential only to determine the maximum 
of the emission barrier. The theory cuts off the image potential short of the sheath and 
plasma (ref. 12). Although the old effective emitter served well, it poses two problems 
in sheath calculations: The absolute overall voltages through the sheath differ for 
opposite charges, and net current results in an isothermal system. These difficulties 
created the need for the virtual Schottky electrode. 
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With a virtual Schottky electrode, both charges undergo the same absolute sheath 
drop, no net charge flows at equilibrium, and the electrode emits all particles under the 
influence of one effective work function. As figure 3 shows, this effective work function 
is the usual work function altered by the Schottky depression (e[Eg |)^^. The Schottky 
depression lowers the work function in the positive -ion sheath and raises the work function 
in the electron sheath. This effect increases electron emission in a positive-ion sheath 
and ion emission in an electron sheath. Conversely, the effective work function decreases 
ion emission in a positive -ion sheath and electron emission in an electron sheath at the 
virtual Schottky electrode. The latter property departs from tradition. But the amount 
of retarded charge from the virtual Schottky electrode that escapes through the sheath 
into the plasma is the same as for the traditional thermionic emitter. This result stems 
from the new overall sheath drop, which equals the traditional sheath drop minus the 
Schottky depression. 

Thus, the overall emissions of electrons and ions from the virtual Schottky electrode 
are identical with those from the traditional thermionic emitter. Also coincident with 
tradition, the image potential does not extend into the sheath and plasma with the virtual 
Schottky electrode. To repeat, however, the new model gives one absolute sheath drop, 
no net current in isothermal systems, and one effective work function for electron, ion, 
and atom emissions. 

The virtual Schottky electrode handles particle escape from the electrode effectively, 
but what of the particles leaving the plasma? Initially, the new overall sheath drop may 
appear troublesome, but using the old sheath drop for both charges from the plasma 
produces net currents for equilibria. Although the old drop accelerates the right amount 
of one charge to the electrode it returns too many of the retarded particles back to the 
plasma. Correcting this error requires two work functions and two sheath drops. The 
virtual Schottky electrode, though, uses one effective work function and one overall 
sheath drop to yield no isothermal flow of net charge. 

The new sheath drop may also seem too small for the accelerated plasma particle. 
However, the new and the old models deliver identical accelerated current densities to 
the electrode. Furthermore, the new sheath drop is correct for the retarded flow from 
the plasma, and according to earlier arguments, the new theory handles all emissions 
well. Thus, the virtual Schottky electrode provides a sheath model that treats all charge 
flows correctly and conveniently. 


Sheath Current Densities 


The fully expanded equations for sheath current densities appear in reference 1. 
Some of those expressions look a bit alien, but they all came from the equations (2a) 
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to (7a), with additions of effective work functions and exponential cutoffs by sheath 
barriers . 

As figure 3 indicates, cp' in equations (2a) to (7a) is w - (eE-^)^^ for positive-ion 
sheaths and <p + (-eEg) ' for electron sheaths. Equations (2a) to (7a) also change for 
intervening sheaths: 


Positive -ion sheath (AVg is positive): 

( 2b ) j e0 E = W * j eIE = W exp 
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Electron sheath (AVg is negative): 
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With appropriate substitutions from this paragraph, equations (2a) to (7a) become those 
of reference 1. 

The net negative current densities into the plasma also result from the equations of 
the preceding paragraph: 


Positive -ion sheath (AVg is positive): 

J = ^eOE + ^iOP ~ ^iOE exp 


e AV C 


kT e 
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Electron sheath (AVg is negative): 


/e AV C 


J = 


jeOE ex P 


kT 


+ j i0P exp 


E 


/e AV c 


* T iP, 


I - j iOE " ^eOP 


(10) 


8 



Of course, these current densities prevail through the respective sheaths. But the 
expressions for current densities within the sheaths are more complicated. Because 
such equations serve only for computing charge densities in the sheaths, they appear only 
implicitly in the next section. 


Sheath Charge Densities 


If the separate charge -flow terms at AV within f he sheath are divided by their 
average velocities, the resulting sum is the charge density at AV. Thus, the charge 
density relations evolve from equations (9) and (10): 

Positive -ion sheath (AV is positive): 


The first term j e Q E in equation (9) represents the flow of electrons from the 
electrode accelerated through AV g - AV to produce an average velocity 
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The second term i i0 p ^ equation (9) gives the flow of ions from the plasma ac- 
celerated through AV. By analogy with equations (11) and (12), the net negative charge 
density becomes 
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The third term j^Qg exp (-e AVg/xTg) in equation (9) results from the flow of elec- 
trode ions through the retarding field of the sheath. At AV, this retardation amounts to 
AVg - AV. Thus, a cutoff flow j-Qg exp £-e(AV g - AV)//cTgj passes AV climbing 

the sheath barrier with an average velocity (2/cTg/ffrm)'^. This part contributes a 
negative charge density: 
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and also effects the negative charge density 


The overall influence of 
equations (14) and (15): 
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The fourth term j e Qp exp (-e AVg//cT e p) in equation (9) indicates the retarding effect 
of the positive -ion sheath on plasma electrons. By inspection of equation (16), the effect 
of this retarded electron stream on the negative charge density at AV is 
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Adding equations (12), (13), (16), and (17) yields the net negative charge density at AV: 
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When AV equals zero, the third term of equation (18) approximates half of the 
positive charge density in plasma. This simple relation allows an initial estimate of 
AVg. 
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Electron sheath (AV is negative): 


Using figure 3, one can analogize the expression for net negative charge density in 
an electron sheath from equation (18): 
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Here, at zero AV, the first term approximates half of the electron charge density in the 
plasma and allows an estimate of AVg. 


Sheath Structures 

With the net charge densities given by equations (18) and (19), simple expressions 
yield electrostatic fields and corresponding distances for assigned AV's through the 
sheaths 
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where x = 0 at AV = AVg and x = Xp = Xg at AV = 0. Although equations (20) and (21) 
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appear simple, their evaluations require numeric integrations. These results are sums 
of the areas of trapezoids based on 20 equal increments of AV from zero to AVg. 
Equation (20) poses no problem; for the machine, neither does equation (21). But at 
AV = 0, E is zero, and equation (21) lacks definition. The machine, however, merely 
uses half of the E at AV/20 as the average for the first increment and computes a 
finite Xg. This distance is an effective sheath width, which includes most of the sheath 
structure. The value of the effective sheath width derives from the extent that it repre- 
sents a fully developed sheath. As reference 1 shows, for near -equilibria and small 
net transport, the positive and negative charges of the sheath almost balance at Xg. 
Thus, the effective sheath width contains practically all the sheath functions within a few 
emission Debye lengths (ref. 1). 


Computing Method 

The iterative sheath solution follows an almost classic numeric method: 

(1) Run screening tests to select noncondensing plasmas that produce essentially 
collisionless sheaths. These tests require effective plasma pressures lower than the 
vapor pressure of the plasma chemical at the electrode temperature and minimum mean 
free paths longer than the Debye length. Reference 1 discusses this screening for 
cesium. 

(2) Solve for current densities deleting effects of sheath fields (eqs. (2a) to (7a)). 

(3) Estimate the overall sheath potential drop AVg (see section Sheath Charge 
Densities). 

(4) Divide AVg into 20 increments and compute charge densities for increasing 
potentials through the sheath. 

(6) Integrate the charge density function over AV to obtain the sheath field (eq. (20)). 

(6) Use the sheath potential drop and field at the virtual Schottky electrode to calcu- 
late new emissions. 

(7) Repeat the cycle until no boundary current density changes by more than 0. 1 per- 
cent of the smallest boundary current density. 

(8) Integrate the negative reciprocal of sheath field over voltage to determine sheath 
distances (eq. (21)). 

(9) Compute and record incremental and overall sheath results for the final iteration. 
The details of these equations, computing methods, FORTRAN statements, and machine 
outputs appear in reference 15. 


linn 


ii n 1 1 ■ i 1 1 it it • hi 


DISCUSSION OF RESULTS 

This report aims primarily at approximating some possible internal potential struc- 
tures for cesium diodes. But on the way to that goal, the results also allow correlations 
of overall sheath characteristics. These generalizations make up figures 4 to 8. Motive 
diagrams and equations in the previous section amplify the meanings of the parameters 
and properties of the sheaths. Because all but the ordinates of figures 4 to 8 come 
directly from the assigned parameters (ref. 1), the reader can estimate sheath conditions 
with these plots rather than with the complicated iterative solutions. 

In figures 4 to 8, the data correlate well because of the selection of the variables: 
Reference 1 and the previous section define the effective sheath width Xg. The emission 
Debye length X^ arose because emission energy rather than plasma electron temper- 
ature dominates charge separation in these sheaths (ref. 1). For this reason, the sheath 
properties depend on emitter temperature Tg and plasma electron concentration N g p. 
The hybrid emission Debye length is a function of these basic variables 

Xj-)g « 6. 9 (Tg/N e p) / (22) 

This emission Debye length correlates the sheath data better than the ordinary Debye 
length. The ratio of the effective sheath width x g to the emission Debye length Xp E is 
the dimensionless sheath width of figure 4. 

In figures 4 to 8, the abscissas are a dimensionless variable termed the Richardson - 
Dushman sheath drop or R-D drop. This R-D drop results from dividing the work function 
minus the plasma potential by the voltage equivalent of the emitter temperature. The 
ordinate of figure 5 is the actual sheath drop divided by the dimensional R-D drop 
(cp - (p ). The dimensionless electrode field of figure 6 is the product of the actual 
electrode field and the emission Debye length divided by the dimensional R-D drop. And 
finally, figures 7 and 8 present the majority charge density at the electrode and the cur- 
rent density for these sheaths. 

In reference 1, two conditions restrict the applicability of the sheath model: Net 
current densities remain below 6 percent of the random particle fluxes in the plasmas. 

And minority charge densities at the effective sheath edge approximate their values in 
the neutral plasmas to within 1 percent. In the present report, however, wider regions 
of sheath similarities result from using less demanding limits: Net current densities 
range up to 11 percent of the random plasma circulations. And minority charge densities 
at the sheath, plasma interface approach their plasma values to within 2 percent. 

This shift constitutes no change in the interpretation of the restrictions of the 
theory - the smaller the net transport and the sheath -edge charge discrepancy, the 
better. The expansion of the bounds merely extends the results to systems of higher 
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diode-power production, hence greater interest. In most of these extrapolations, how- 
ever, the sheaths yet resemble equilibrium and near -equilibrium versions from refer- 
ence 1, as figures 4 to 8 reveal. 

Effective sheath widths for separated electrode and plasma temperatures are shown 
in figure 4 to be close to the equilibrium correlation at high R-D drops. While the re- 
sults for electrode temperatures lower than those of the plasmas behave like those of 
reference 1, the curve for electrode temperatures higher than those of the plasmas 
breaks upward as R-D drops decrease. This departure occurs because the net electron 
current grows exponentially with electrode temperature. Thus, the minority charge 
strongly influences the structure of these positive -ion sheaths at low overall potential 
drops. When the electron current imbalance comes from raising the plasma temperature, 
though, the outcome approximates that for increasing the electron temperature only at 
the same electron concentration (refs. 1 and 5). 

In figure 5, the present overall sheath drops appear quite similar to those for 
equilibria. The emitter fields of figure 6 duplicate the equilibrium trends and extend 
them into straight semilogarithmic lines. These asymptotic exponentials are approxi- 
mately three times the R-D drop. As figure 7 illustrates, the majority-charge density 
at the emitter again depends only on the plasma electron number and the R-D drop. And 
the net current densities given in figure 8 look much like those of reference 1. 

In general, sheaths between an electrode and a cesium plasma with either an 
elevated or a reduced temperature compare with those for equilibria and for systems with 
only the electron temperatures increased (ref. 1). Because changing plasma and plasma- 
electron temperatures produces similar effects, decreasing the plasma -electron temper- 
ature should also yield results like those of figures 4 to 8 and of reference 1. 

Returning from generalization to the specific goal, figure 9 presents matched pairs 
of electrode sheaths in cesium plasmas. These coupled sheaths pass net negative current 
densities, which match within 2 percent, from the emitter to the collector. Complete 
tabulations of detailed and overall results attend each sheath; definitions of symbols ap- 
pear in appendix B. In addition to the points for the graphic potential distributions, these 
tables give data for profiles of charged-particle number densities and of electric fields 
in thermionic diodes. 

Although figure 9 shows no changes in potential across the plasmas, separating the 
zero potential levels of corresponding emitter and collector sheaths by a suitable amount 
approximates the plasma drop (refs. 2 to 4). For these relatively small net current 
densities, the necessary plasma fields are low compared with those at the electrodes. 
Such gradual gradients in the plasmas comply with the assumptions of the sheath model. 

Because reference 1 dwells on the limitations of the theory, this work mentions 
them only in passing. But one problem requires emphasis: In practice the plasma, 
sheath, and electrode determine the work function; it is not an assignable parameter. 
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Thus, though the systems given here might be approximated by special selections, they 
serve as examples, not recommended diodes. Yet the demonstrated generalizations of 
sheath effects make the specific results practical as well as heuristic. To reveal proper- 
ties for particular electrode materials, however, calculations must include variations of 
work functions with cesium arrival rates. 

Lacking this refinement and its computational complication, the present studies still 

yield emission barriers not far removed from reality. The work functions vary for 

emitters from 3. 8 to 4. 9 volts and for collectors from 3. 5 to 4. 1 with ion sheaths and 

13 

from 2.0 to 2. 2 with electron sheaths. In the cesium plasmas with 10 electrons per 

cubic centimeter, the net current densities range from 0. 094 to 1. 2 amperes per square 

12 

centimeter through effective sheath widths from 1. 9 to 2. 3 microns. At 10 electrons 
per cubic centimeter, 0.010 to 0.095 amperes per square centimeter pass through 
5. 8- to 6. 9-micron sheaths. 

Such widths take up little of the typical 250-micron gap between electrodes in a 
thermionic converter. But the fields and potential changes of the sheaths strongly in- 
fluence diode performances. Overall sheath drops in this work range from 0. 83 to 
1. 8 volts for emitters and from 0. 35 to 0. 82 in absolute values for collectors. And if 
small plasma losses occur, these cesium cells produce approximately 0. 01 to 1 watt 
per square centimeter. 

Although interesting diodes usually generate more than 10 watts per square centi- 
meter, the lower levels given here represent a part of the spectrum of productive energy 
conversion. Furthermore, the conditions depart considerably from those for the 
isothermal diode (refs. 1 and 12). The data, therefore, extend the theoretic region 
available for comparison with experiment and prediction of performance in thermionics. 

In addition, the generalizations of overall characteristics for the sheath model of refer- 
ence 1 should add to the mechanistic understanding of plasma energy converters. 


CONCLUDING REMARKS 

These low -current approximations of cesium diodes dramatize sheath effects. Al- 
though the interelectrode space is small enough to be a mechanical problem, the effective 
emitter and collector sheaths take up little of the diode gap. Yet the sheath drops are 
similar to diode output voltages. Little sheath widths and large potential changes mean 
big fields at the electrodes. These fields depress emission barriers for particles ac- 
celerated from the electrodes and swell current densities. Because diode emitters often 
develop positive -ion sheaths, their fields produce electron emissions greater than basic 
work functions indicate. And because ion sheaths increase the cesium delivery to the 
electrode, the basic work functions run lower than those indicated by cesium arrival 
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rates based on plasma circulation. The lower work functions also yield more electron 
emission. Thus, positive-ion sheaths contribute to greater diode currents with two 
characteristics: intense fields and high cesium arrival at the emitter. 

Although transport within the sheath affects diode outputs directly, these sheath 
effects on electron emission appear more important. Many research programs now seek 
emitters with low work functions at high temperatures to boost emission. This problem 
is difficult, particularly because low work functions often mean high vapor pressures and 
added impurities. In the face of this need and its cost, the idea of improving emission 
with a strong positive-ion sheath looks promising. The ion sheath uses the electrodes 
and the plasma that are the basic diode to increase fields and cesium concentrations at 
the emitter. When both these sheath properties lower the effective work function of 
the emitter, electron emission and diode power grow rapidly. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 21, 1967, 

120-33-02-01-22. 
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APPENDIX A 


SYMBOLS 


E electric field 

e electronic unit charge 

h Planck's constant 

I ionization potential 

J overall net current 

j current density or net particle 

current density 

m particle mass 

N particle number density 

T temperature 

V potential 

AV potential relative to plasma potential 

v velocity 

x distance from emitter 

k Boltzmann constant 

fx Fermi level 

emission Debye length, 

(f work function 

<£•' effective work function for virtual 

Schottky emitter 

(p work function (plasma potential) 

that precludes an ion sheath 
(used for plasma potentials for 
both ion and electron sheaths 
where they appear on the same 
figure) 


(p QO work function (plasma potential) 
that precludes an electron 
sheath 

Subscripts: 

a atom 

e electron 

E emitter 

i ion 

IE into emitter 

IP into plasma 

O overall (from work function to 

plasma potential) 

OE out of emitter 

OP out of plasma 

P plasma 

S across sheath 

AV at AV 
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APPENDIX B 


FORTRAN 

symbol 

I 

TE 

PHI 

NEP 

TEP 

TIP 

LAMBDA 

PV 


NOMENCLATURE FOR IBM OUTPUT 

Description 


Algebraic 

symbol 


I 

T 

<P 


E 


N eP = N iP 


eP 


T iP ” T aP for these 


solutions 


DT 


6 . 


eP' 


\l/2 


eP 


VN 


eP/ 


vp 


ionization potential for plasma atoms 
emitter temperature 
work function 

plasma electron number density 
plasma electron temperature 
plasma ion temperature 

plasma Debye length 

vapor pressure of plasma element at 

T-c 


-3 


Units 

V 
°K 

V 
cm 
°K 
°K 

cm 


torr 


N/m 




1/2 ‘ \ 

1.333x10 

LAMBDA(TE) 

X »6.9(-E- 

E \ N eP, 

•) emission Debye length 

cm 

DV 

AV 

sheath potential measured from plasma 
electron potential 

V 

ND(DV) 

N e AV “ N i AV 

net number density of charged particles 
at AV 

-3 

cm 

NE(DV) 

N e AV 

electron number density at AV 

cm - ^ 

NI(DV) 

N i AV 

ion number density at AV 

cm"^ 

E(DV) 

E AV 

electron electrostatic field at AV 

V/cm 

X(DV) 

X AV 

distance from emitter to AV 

cm 

JEE 

•'eE 

emitted electron current density 

A/cm 2 

JEP 

■*eP 

plasma electron random current 
density 

A/cm 2 

JIP 

jip 

plasma ion random current 
density 

A/cm 2 
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FORTRAN Algebraic 

symbol symbol 


JAP 

■*aP 

J 

J 

PP 

Pp 

JIE 

•*iE 

JAE 

^aP 

JA 

•*a 

JI 

h 

JE 

j e 

JA/JAP 

^a^ ^aP 

JE/JEP 

•^e^eP 

JI/JIP 

j i/j iP 

DVS 

AV S 

XDVS 

X AV S = X P= X J 

NAP 

N aP 

XD/LAM 

x s A d 

SC 

<- E ) 1/2 

PHZ 

Vo 

EDVS 

E av s - E E 

DVSRD 

0 

1 

11 

o 

> 

< 

DVS/RD 

AV s _ AV s 
AV o * - % 

ELM/RD 

e e x d // ^ " 


Description 

Units 

plasma atom equivalent random 

A/cm 2 

current density 


net current density through sheath 

A/cm 2 

plasma pressure 

torr 

/ o' 


N/m 2 


1.333X10 2 


2 

emitted ion current density A/cm 

2 

emitted equivalent atom current A/cm 

density 

2 

net equivalent atom current density A/cm 

2 

net ion current density A/cm 

2 

net electron current density A/cm 


overall sheath voltage drop V 

effective sheath thickness cm 

plasma atom number density cm 


Schottky depression of work function V 

plasma potential (work function for V 

no sheath) 

electrostatic field at emitter V/cm 

Richardson-Dushman overall sheath V 

voltage drop 
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FORTRAN 

PHZZ 

DRD/KT 

NTP 

NCE 

NTE 

RD/KTE 

X/LMTE 

ELT/RD 

NEPA 

NIPA 


Algebraic 

symbol 

c P (<p . = a, for 
oo oo ^o 

equilibrium and 
electron sheaths) 

e W ~ <P Q l/* T e 


Description 


Units 


e k - <p 0 IA t e 

X s/ X DE 

e e x de // ^ " < P 0 ) 


plasma potential at equilibrium 

(work function for no sheath and no 
net current) 

-3 

total particle number density in plasma cm 

-3 

total charge number density at emitter cm 

-3 

total particle number density at emitter cm 

-3 

plasma electron number density from cm 

sheath calculations 

3 

plasma ion number density from sheath cm 
calculations 
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(a) Positive-ion sheath* For electrons, E and AV are positive; ef- 
fective work function, <p' =<p- (-eE E )^ 2 . 


Electrode Motive Plasma 



(b) Electron sheath. For electrons, E and AV are negative; effec- 
tive work function, <P‘ = <P + (-eE E )^ 2 . 


Figure 3. - Plane sheaths between virtual Schottky electrodes and near- 
equilibrium plasmas. 
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Figure 4. - Effective sheath widths for emitters and collectors in cesium plasmas. 
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Figure 7. - Charge concentrations at emitters and collectors in cesium plasmas. 
At |e(#> - (p 0 )/k.T e | * 0, N eE - N ep ■ N iE - N ip . 
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Net negative current density, A/cm 



Figure 8. - Charge transport through emitter and collector sheaths in cesium 
plasmas. 
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(a~l) Current density, 1. 16 amperes per square centimeter. 

(a) Emitter at 2400° K, plasma at 2000° K with 10 13 electrons per cubic centimeter, and collector at 1600° K. 
Figure 9. - Some sheath structures for cesium thermionic diodes. 
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(a-1) Concluded. Current density, 1. 16 amperes per square centimeter. 

(a) Continued. Emitter at 2400° K, plasma at 2000° K with 10 13 electrons per cubic centimeter, and collector at 1600° K. 

Figure 9. - Continued. 
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(a-2) Current density, 0. 565 ampere per square centimeter. 

(a) Continued. Emitter at 2400° K, plasma at 2000° K with 10 13 electrons per cuoic 
centimeter, and collector at 1600° K. 


Figure 9. - Continued. 
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(a-2) Concluded. Current density, 0. 565 ampere per square centimeter. 

(a) Continued. Emitter at 2400° K, plasma at 2000° K with 10 13 electrons per cubic centimeter, and collector at 1600° K. 

Figure 1 ?. - Continued. 
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(a-3) Concluded. Current density, 0.0951 ampere per square centimeter. 

(a) Concluded. Emitter at 2400° K, plasma at 2000° K with 10 13 electrons per cubic centimeter, and collector at 1600° K. 

Figure 9. - Continued. 
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(b-1) Concluded. Current density, 1. 17 amperes per square centimeter. 

(b) Continued. Emitter at 2200° K, plasma at 2000° K with 10 13 electrons per cubic centimeter, and collector at 1800° K. 

Figure 9. - Continued. 
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(b-2) Current density, 0. 568 ampere per square centimeter. 

(b) Continued. Emitter at 2200° K, plasma at 2000° K with 10 13 electrons per cubic centimeter, and collector at 1800° K. 

Figured - Continued. 
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(b-2) Concluded. Current density, 0. 568 ampere per square centimeter. 

(b) Continued. Emitter at 2200° K, plasma at 2000° K with 10 13 electrons per cubic centimeter, and collector at 1800° K. 

Figure 9. - Continued. 
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<b~3) Current density, 0, 0945 ampere per square centimeter. 

(fa) Continued. Emitter at 2200° K, piasma at 2000° K with lO 1 ^ electrons per cubic 
centimeter, and collector at 1800° K. 

Figure 9. - Continued. 
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(b-3) Concluded. Current density, 0. 0945 ampere per square centimeter. 

(b) Concluded. Emitter at 2200° K, plasma at 2000° K with 10 i3 electrons per cubic centimeter, and collector at 1800° K. 

Figure 9. - Continued. 
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(c-1) Current density, 0. 0939 ampere per square centimeter. 

(c) Emitter at 20(XT K, plasma at 1800 n K with electrons per cubic centimeter, and collector at 1600° K. 


Figure 9. - Continued. 
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(c-1) Concluded. Current density, 0. 0939 ampere per square centimeter. 

(c) Continued. Emitter at 2000° K, plasma at 1800° K with 10 12 electrons per cubic centimeter, and collector at 1600° K. 

Figure 9. - Continued. 
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ic-2) Current density, 0. 0601 ampere per square centimeter. 

(c) Continued. Emitter at 2000° K, plasma at 1800° K with 10 12 electrons per cubic centimeter, and collector at 1600° K. 

Figure 9. - Continued. 
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(c-2) Concluded. Current density, 0. 0601 ampere per square centimeter. 

(c) Continued. Emitter at 2000° K, plasma at 1800° K with 10 12 electrons per cubic centimeter, and collector at 1600° ft 

Figure 9. - Continued. 
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(c-3) Current density, 0. 00994 ampere per square centimeter. 

(c) Continued. Emitter at 2000° K, plasma at 1800° K with 10 12 electrons per cubic centimeter, and collector at 1600° K. 


Figure 9. - Continued. 
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(c-3) Concluded. Current density, 0. 00994 ampere per square centimeter. 

(c) Concluded. Emitter at 2000° K, plasma at 1800° K with 10 12 electrons per cubic centimeter, and collector at 1600° K. 

Figure 9. - Continued. 
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(d-1) Current density, 0. 0945 ampere per square centimeter. 

(d) Emitter at 2000° K, plasma at 1700° K with 10 12 electrons per cubic centimeter, and collector at 1400° K. 

Figure 9. - Continued. 
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(d-1) Concluded. Current density, 0. 0945 ampere per square centimeter. 
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(d) Continued. 


Emitter at 2000° K, plasma at 1700° K with 10 12 electrons per cubic centimeter, and collector at 1400° K. 


Figured - Continued. 
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(d) Continued. Emitter at 2000° K, plasma at 1700° K with 10 12 electrons per cubic 
centimeter, and collector at 1400° K. 
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(d-2) Concluded. Current density, 0. 0603 ampere per square centimeter. 

(d) Continued. Emitter at 2000° K, plasma at 1700° K with 10 12 electrons per cubic centimeter, and collector at 1400° K. 

Figured - Continued. 
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(d-3) Current density, 0. 0101 ampere per square centimeter. 

(d) Continued. Emitter at 2000° K, plasma at 1700° K with 10^ electrons 
per cubic centimeter, and collector at 1400° K. 


Figure 9. - Continued. 
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(d-3) Concluded. Current density, 0. 0101 ampere per square centimeter. 

(d) Concluded. Emitter at 2000 ? K, plasma at 1700° K with 10 12 electrons per cubic centimeter, and collector at 1400° K. 

Figure 9. - Concluded. 
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